Effects of the interplay between interaction and disorder in bipartite entanglement 
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We use a disordered anti-ferromagnetic spin-1/2 chain with anisotropic exchange coupling to 
model an array of interacting qubits. All qubits have the same level spacing, except two, which are 
called the defects of the chain. The level spacings of the defects are equal and much larger than all 
the others. We investigate how the entanglement between the two defects depends on the anisotropy 
of the system. When the anisotropy coupling is much larger than the energy difference between a 
defect and an ordinary qubit, the two defects become strongly entangled. Small anisotropics, on the 
contrary, may decrease the entanglement, which is, in this case, also much affected by the number 
of excitations. The analysis is made for nearest neighbor and next-nearest neighbor defects. The 
decrease in the entanglement for nearest neighbor defects is not very significant, especially in large 
chains. 
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I. INTRODUCTION 

Entanglement is one of the most striking properties of 
quantum mechanics. It describes a correlation between 
quantum mechanical systems that does not occur in clas- 
sical mechanics, namely that a pure state of a composite 
quantum system cannot be written as a product of the 
states of its constituents 0. It has long been a subject 
of interest in foundations of quantum mechanics, but re- 
cently it has acquired a fundamental role in quantum 
computation and quantum information . 

Several attempts to quantify entanglement have been 
developed 0, but this is an area of research still in 
progress. In the case of a pure state of a bipartite system, 
a good and widely accepted measure of entanglement is 
the von Neumann entropy of its reduced density matrix. 
This entropy can be associated with a quantity called 
concurrence which we adopt here as our measure of 
entanglement. The concurrence varies from to 1, the 
larger it is the more entangled is the state. 

Spin chains are ideal systems for the study of entan- 
glement. They are naturally used to model quantum 
computers (QC's): the two states of a spin-1/2 particle 
correspond to the two levels of a qubit and the exchange 
interaction corresponds to the qubit-qubit interaction. A 
particle with spin up corresponds to an excitation or an 
excited qubit. We consider a spin chain described by the 
XXZ model. This is the model used in quantum com- 
puters based on electrons on helium The Ising part 
of the interaction is proportional to the anisotropy cou- 
pling, as shown in the next section. It only plays a role 
when two or more excitations are present in the chain. 
The XY part of the interaction is responsible for hopping 
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the excitations and can be used to create entanglement 
between two or more qubits |^, Q ■ However it can also 
quickly destroy bipartite entanglement, because the ex- 
citations soon spread all over the chain. 

A way to avoid the spreading (or delocalization) of the 
excitations is by tuning the qubits away from resonance 
■ Such control of the qubit level spacings enables us to 
entangle just some specific qubits. This is the strategy 
we use in this paper. To have two chosen qubits maxi- 
mally entangled, we tune them in resonance. Their level 
spacings are different from all other qubits, they are the 
defects of the chain. 

When there is only one excitation in the chain, the en- 
tanglement between the two selected sites becomes triv- 
ial. By assuming that the energy difference between the 
defects and the ordinary qubits is much larger than the 
strength of the XY interaction, the Hamiltonian of the 
system gives two maximally entangled states. They cor- 
respond to superpositions of the basis states where the 
excitation occupies one of the defects 0- However, when 
more than one excitation is present, the entanglement 
between the defects can be significantly affected by the 
Ising interaction. The main purpose of the present work 
is to analyze which is, in this case, the eigenstate with 
maximum concurrence and how the value of the maxi- 
mum concurrence depends on the anisotropy. 

We find numerically that for large values of the ani- 
sotropy coupling, nearest neighbor defects are strongly 
entangled and the value of the concurrence is almost neg- 
ligibly affected by the number of excitations in the sys- 
tem. In the case of separated defects, this observation is 
valid only for long chains. In the opposite situation of 
small anisotropics, the concurrence between the defects 
is strongly dependent on the number of excitations and it 
can become quite small. Partial analytical justifications 
for these observations are provided. 

The paper is organized as follows. In Sec. II we de- 
scribe the model and the adopted measure for bipartite 
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entanglement. In Sec. Ill we analyze the maximum 
concurrence between nearest neighbor and next nearest 
neighbor defects for various values of the anisotropy cou- 
pling. A discussion on how to prepare maximally en- 
tangled states and how the anisotropy may help is also 
presented in this section. 



II. THE MODEL AND THE MEASUREMENT 
OF ENTANGLEMENT 

We consider a spin chain with nearest neighbor inter- 
action. The Hamiltonian describing the system is given 
by the XX Z model 



hop 



n=l 



(1) 



where fi = 1 and a^-'^'" are Pauli matrices. There are 
L sites and we deal with a periodic (or closed) chain, 
that is, sites n + L and n are the same. Each site n 
is subjected to a magnetic field in the z direction, giving 
the energy splitting e„. In this description, the excitation 
energy of a qubit is the Zeeman energy of a spin. In this 
disordered chain, not all the qubits have the same level 
spacing e. The energy of two qubits, called defects, differs 
from those of other qubits by d. The parameter J is 
the hopping integral and A is a dimcnsionless parameter 
related to the anisotropy coupling. The diagonal term 
a^a^_^_i gives the Ising interaction and the non-diagonal 
term Tfhop is responsible for propagating the excitations. 
We set d, J and A > 0. 

The anisotropy in the Hamiltonian Q is different from 
the one considered in some previous models |^ ITol llll | . 
The hopping part of the Hamiltonian can be equiva- 
lently written as i/hop oc J ^u'^^u'^^j^^ -I- Jy^n'^^^x- Here 
Jx = Jy = J, but in the models cited above, the degree 
of anisotropy comes from the difference between and 
Jy. The anisotropy in our case originates from the extra 
Ising interaction. In terms of entanglement, few studies 
have been developed with this model 0, 1^. More- 
over, our goal is not simply to analyze how entanglement 
may depend on the described anisotropy, but we aim to 
analyze how it is affected by the interplay between this 
interaction and disorder. 

The disorder of the system we consider is characterized 
by the presence of the two defects. In principle, total con- 
trol of the qubit level spacings is available, which allows 
the creation of defects. 

The Ising part of the XX Z model combined with de- 
fects has been used before to create maximally entangled 
states 0. The study of entanglement with impurities, 
but without any (T^cr^+i interaction, has also been done 
in [lol [Tll |. We emphasize the importance of studying 
the effects of this extra interaction in disordered systems, 
which are far from trivial. The Ising interaction is only 



relevant when at least two excitations are present, be- 
ing therefore associated with many-body problems. It is 
actually at the heart of one of the most challenging prob- 
lems in condensed matter physics, namely the difficulty 
in localizing many-particle states 

In what follows, we count energy off the ground state 
energy Eq = —{Ls + 2d)/2 + LJA/4, i.e., we replace in 
Eq. ^ H H - Eq. 

To address the different states of the system we use a 
notation that is common in the study of spin chains with 
the Bethe ansatz The state corresponding to one 

single excitation on site n, that is | J,i|2 ■■• in-iTnin+i 
... Il), or equivalently |0i02...0„_i1„0„+i...0l), is simply 
written as (j>{n). The state of two excitations, one on 
site n and the other one on site m, is (j){n,m), which is 
a simplified notation for | I1I2 ■•• Tnin+i ■■■ Tm ••■ II 
), or equivalently |0i02...f nOri+i...lm---OL)- Basis states 
where each excitation is confined to a single site, such 
as (/)(n,m), are called in quantum computing quantum 
registers. These are the states where the measurements 
are performed. 

Since we want to study the effects of the Ising interac- 
tion, several excitations have to be considered, which lim- 
its the numerical analysis to small chains. In the model 
described by Eq. the z component of the total spin 
is conserved, so states with different number of 
excitations are not coupled and the Hamiltonian is made 
of uncoupled blocks. However the blocks can still be very 
large. For a chain with L sites, the block corresponding 
to N excitations has dimension Ll/[Nl{L — N)l]. Its di- 
agonalization leads to eigenstates which correspond to 
linear superpositions of quantum registers (j) with N ex- 
citations. Each state k is written as 



\^nI)= E a'^''Hn,m,...,Nmn,m,...,N)). 

n<m<...<N=l 

(2) 

In order to study quantitatively the entanglement be- 
tween two qubits we calculate their entanglement of for- 
mation Ep. Given the density matrix p that describes 
our pair of qubits, Ep is the average entanglement of the 
pure states of the decomposition of p, minimized over all 
possible decompositions: 



Ef{p) = min^ piE{ipi), 



(3) 



where J2iPi = 1, < < 1, and p = J2iPi IV"*) (V"*!- 
Here E{ip) is the von Neumann entropy of either of the 
two qubits 0. Wootters et al 3 have shown that, for 
a pair of qubits, Ep is a monotonically increasing func- 
tion of the concurrence, which one can prove to be an 
entanglement monotone. Since the concurrence is math- 
ematically simpler to deal with than Ep, we adopt it here 
to measure the entanglement between two qubits. It is 
given by 0: 



C = max{Ai — A2 — A3 — A4, 0}, 



(4) 
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where Ai, A2, A3, and A4 are the square roots of the eigen- 
values, in decreasing order, of the matrix R = pp. The 
matrix p is the time reversed matrix 



/) P* K 



(5) 



The symbol p* means complex conjugation of the matrix 
p in the basis {|11) , |10) , |01) , |00)}. 

Maximum entanglement corresponds to C = 1 and no 
entanglement gives C — 0. To compute the concurrence 
of two qubits in a chain with several sites, we trace over 
the qubits we are not interested in and study the reduced 
density matrix of the two chosen ones. 



III. NUMERICAL RESULTS AND 
ANALYTICAL ANALYSIS 

We study the entanglement between the two defects of 
the chain. The difference in energy between them and the 
other qubits, d, is assumed much larger than the hopping 
integral, d >> J. This guarantees that an excitation 
placed on one defect can only hop between the two defect 
sites. As a consequence, among all the eigenstates of the 
system, the ones with a single excitation shared between 
the defects have the largest concurrences. 

In the first subsection below, we analyze the case where 
the defects correspond to two nearest neighbor qubits and 
in the second subsection they are two next nearest neigh- 
bors. Since the chain has periodic boundary conditions, 
any pair of qubits (no, mo) translated through the chain 
is equivalent, so we choose, for the numerical calcula- 
tions, the pair 1 and 2 as nearest neighbor defects and 
the pair 1 and 3 as next-nearest neighbor defects. 

We study how the maximum concurrence between the 
defects depend on the anisotropy and which is the corre- 
sponding eigenstate. 
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FIG. 1: (color online) Maximum concurrence vs. anisotropy 
coupling in closed Heisenberg spin-1/2 chains with defects lo- 
cated on sites 1 and 2. The energy difference between the 
defects and the other qubits is d = lOJ. The panels corre- 
spond to data for chains of different lengths: L = 8 (top), 
L = 10 (middle) and L = 12 (bottom). The maximum num- 
ber of excitations considered for each chain is L/2. The curves 
are chosen as follows: solid (black) curves for 2 excitations, 
dotted (red) for 3 excitations, dashed (blue) for 4 excitations, 
long dashed (green) for 5 excitations, and dot-dashed (violet) 
for 6 excitations. 



A. Nearest neighbor defects 

The numerical results for the dependence of the maxi- 
mum concurrence, Cmax, on the anisotropy of a periodic 
Heisenberg chain with two neighbor defects are shown in 
Fig.n When A = 0, all states with one excitation on the 
defects have Cmax — 1- For small A's, in general, Cmax 
decreases with the number of excitations. The compar- 
ison of Cmax for chains of different sizes but with the 
same number of excitations indicates that, in most cases, 
smaller chains are more affected by the Ising interaction. 
When there are just two excitations, the minimum value 
of Cmax happens when JA = d, while in the case of more 
excitations, this occurs for smaller values of A. When- 
ever JA 3> d, the maximum concurrence stabilizes in a 
value close to 1 (the larger the chain the closer to 1 it 
will be). In this case, the dependence on the number of 
excitations becomes little noticeable. In what follows, we 
try to find justifications, some times analytically, to these 
observations. 



The eigenvalues and eigenstates of an anisotropic spin 
chain with no defects can be analytically obtained with 
the Bethe ansatz In the presence of defects, where 
d ~ J, most commonly the chain becomes non integrable 
[T^ Ir^ . When d » J, as in this paper, we are again 
capable of solving the eigenvalue problem in certain sit- 
uations. This happens because, with such large d, the 
XY-type interaction between defect and ordinary qubit 
is negligible, there is no hopping of excitations between 
them. It works as if we had cut the closed chain and cre- 
ated two chains with free boundaries (two open chains): 
a small one corresponding simply to the two defects and 
another one with L — 2 sites. 

In the simple case of a single excitation in the chain, 
the spectrum is divided into two well separated bands. 
One is made of L — 2 states with energy in the range 
El — J < El < El + J, where ei — e — J A. They are 
states with the excitation out of the defects. The other 
band is made of just the following two linear superposi- 
tions [0(no) ± (/)(no -I- l)]/^/2, with energy ei -|- d ± J/2, 
respectively. They correspond to Bell (or EPR) states 
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and have concurrence equal to 1. 

In this paper we want to analyze how the Ising inter- 
action affects the entanglement and so we consider two 
or more excitations. 



1. Two excitations and A = 

In the case of two excitations, but with A = 0, we 
should obtain results similar to ones for one excitation, 
that is the states with maximum concurrence should give 
Cmax = 1, since no Ising interaction exists. The states of 
maximum concurrence between the defect sites have one 
excitation shared between them. They are superpositions 
of the quantum registers with energy £2 + d, where £2 = 
2£ — 2 JA and form an energy band that we will call defect 
band, or d-band for short. The transitions between these 
quantum registers are schematically shown below. Each 
quantum register (apart from the border ones) is coupled 
to three others: two corresponding to an excitation hop 
to an ordinary qubit and one associated to an excitation 
hop to the other defect. 



no, no + 2 
I 

no + l, no + 2 



no, no + 3 

I 

no + l, no+3 



no,no + L- 1 

I 

... no + l,no+L-l 



The states of the d-band have the form 



IV': 



(fel,fc2)\ 

2,L / 



= ^c{n,m)\(p{n,m)) 



no — 1+i 



^ c{no,m)\(j){nQ,m)) 

m— no+2 
no-l+L 

+ ^ c{no + l,m)\<p{no + l,m)). (6) 

m— no+2 

[To avoid a heavy notation, the superscript (fci, fc2) was 
suppressed from the right hand side.]. 

We use the Bethe ansatz method to find the coeffi- 
cients c(n, m) for these 2(L — 2) states. The Schrodinger 
equation for c(n, m) is 



E2c{n,m) = {e2 + d)c{n,m) + ^[c{n — l,m) 

+ c{n + 1, to) -|- c(n, TO — 1) + c(n, to + 1)]. 

(7) 

Since the studied chain became equivalent to two open 
chains, it is natural to write 



E2 = e2 + d + Ea + Eb, 
Ea = J cos a, 
Eb = Jcos/3. 



(9) 



Using the boundary conditions for the small chain of de- 
fects 



a(no - 1) = a{no + 2) = 0, 



(10) 



and the boundary conditions for the chain of ordinary 
qubits 

b{na + I) = b{no + L) = 0, (11) 

we obtain 

c(n, to) = Asin[a(n — uq + 1)] sin[/3(m — no — 1)], 

where A is a normalization constant and a and /3 are 
given by 

a = 7rfci/3, fci = l,2; (12) 
(3 = TTk2/{L-l), = l,2,...,i-2. (13) 

Each one of the 2(L — 2) states have then a "Bell-type" 
form 



riQ — l+L 

\^P{n,m)'^'''^''^^) = J2 ^sin 

m— nn+2 



7Tk2{m — riQ — 1) 



X [0(no, to) + (-l)fci+i(/,(,^o + 1, to)] . (14) 

To calculate the concurrence between the two defects, 
we obtain the reduced density matrix of the states above 
by tracing over the ordinary qubits. Only the four ele- 
ments in the middle are different from zero. We have 



/O o\ 

1/2 (-l)'=i+V2 

(-l)'=i+V2 1/2 

\ 0/ 



(15) 



Therefore, when A = 0, all the states in the d-band have 
concurrence equal to 1. These analytical results agree 
very well with the numerical ones. 



2. Two excitations and < A < d/J 



c{n,m) = a{n)h{m) = 

[Aie'"" + ^26"'""] [-Bie"^™ -h 526""^"] . (8) 

The energy E2 as a function of a and [3 is obtained 
from Q, ||HJ) and has the form 



The situation changes as we start increasing A. The 
scheme for the involved transitions from the previous 
sub-subsection still applies, but the quantum registers 
(j){nQ — 1, no) and (^(ng -I- 1, no + 2) now have on-site en- 
ergy equal to £2 -l-d+ JA, which is larger than the energy 
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of the other quantum registers considered. As a conse- 
quence, the coefhcients c{n, m) are not simply written 
as the product 0, which leads to the "Bell-type" states 
given by Ea. ((n|l . 

However, Eq.® can still be a reasonable approxima- 
tion to obtain the state with the lowest energy in the ri- 
band. This state necessarily has small coefficients for the 
exceptional registers (/'(no — 1, uq) and 0(rio + 1, -I- 2), 
so the effects caused by their different energies, which 
move the eigenstates away from the "Bell-type" states, 
are not so significant. Such "ground" state of the band 
has, as expected and also confirmed numerically, the 
largest concurrence. It can be found analytically by us- 
ing Eas.lP^. lfTIUl . and some new boundary conditions for 
the chain of ordinary qubits. This is a relatively good 
approximation to the actual state and provides a good 
description for the decay of Cmax with A. 

Since we are interested in the d-band ground state, 
we select ki = 2, which gives the smallest cos a. The 
boundary conditions for the coefficients b{m) are now 
different. When the first excitation is on site n = no, the 
second excitation, necessarily out of the defects, has to 
satisfy 

6(no + 1) 0. (16) 

Moreover, in this case, the second excitation is also sub- 
jected to the condition 

JAb{7io-l+L) + ^[b{no-2+L)] ^ Ebb{no~l+L), (17) 
which, when combined with the general equation, 

^[bim-l) + b{m + l)]^ Ebbim), (18) 

leads to 



JAb{no-l + L) = -b{nn + L). (19) 

Equivalently, when the first excitation is on site n = 
no -l- 1, the boundary conditions for the second one are 

JAb{no + 2) = ib{no + l), 

b{no + L)^0. (20) 

The wave function in the case of A < d/J is then 
written as 

no~l + L 

\ip{n,m)) = A [S{no + L — m)(j){nQ,m) 

m— no+2 

+S{m-no~l)4>ino + l,m)], (21) 



where A is a normalization constant and the function 
S{x) is given by 

Six) = 2Asm[(3{x - 1)] - sin[(3x]. (22) 

The L — 2 different and non-trivial values of /3 are ob- 
tained from 

2A sin[/?(L - 2)] = sm[(3{L - 1)]. (23) 

We notice once again that we are only interested in the 
value of (3 that gives the smallest cos /?. The eigenvalues 
obtained with Eas.(j^. (|12() . and (|23(l are fairly similar 
to the ones obtained numerically, but the results for the 
eigenvectors are so not good. The approximation 
used here is only acceptable for the lowest-energy state 
of the d-band. Its reduced density matrix becomes 

/O 0\ 

1/2 \A\\s 

^~ \A\^s 1/2 ' ^^^> 

\0 0/ 

where 

na-l + L 

s— S{nQ + L — m)S{m — uq — 1). (25) 

m— no+2 

The maximum concurrence obtained, 



Cmax = ^\A\^S^ + \A\^S+l/i 

-Vl^l^s'- l^ps + 1/4, (26) 

is clearly smaller than 1. \A\'^s can only be equal to 1/2 
when A = and |S'(no + L — m)\ — \S(m — no — 1)|, that 
is when all quantum registers forming the d-band are in 
resonance. However, S'(no + L ~ m) and S{m — no — 
1) differ mostly at the borders of the chain of ordinary 
qubits (where TO = no-l-2orm = no+i— 1) and this 
difference increases as A increases. 

In the case of 8 sites, d/J = 10 and A = 3, for ex- 
ample, we find Cniax ~ 0.91, while numerically we have 
C'max ~ 0.98. For a longer chain with 12 sites and the 
same parameters, Cmax ~ 0.97, while numerically we 
have Cmax ^ 0.99. The agreement between numerics 
and analytics are not excellent, but the approximation 
gives the correct trend: the concurrence decreases as A 
increases (with A < d/J), and the effects become less 
important in larger chains. 

3. Two excitations and A = d/J 

In the case of resonance, when A = d/J, the quan- 
tum registers (f>{nQ + L — l,no) and (/)(no -I- l,no -|- 2) 
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have very large energy and do not participate in the for- 
mation of the eigenstates in the defect band anymore. 
However, quantum registers corresponding to the bound 
pairs (/)(no ~ 2, no — 1) and (f>{nQ + 2, no + 3) now play an 
important role, for they are in resonance with the other 
registers and couple with them in first order of pertur- 
bation theory in J. As a result, a transition from the 
sites (no -I- 1, no + 3) to (no -t- 2, no -I- 3) [equivalently from 
(no — 2, no) to (no — 2, no — 1)], which removes the exci- 
tation from the defect, can happen. This is the cause of 
the drop in the value of the maximum concurrence. 

The scattering process is schematically shown below. 
There are four states out of the core that could lead to a 
Bell-type state. 

no — 1, no + L — 2 
I 

"0i"0 + 2 ^no,no+3^ ... no,nQ + L — 2 

I I 
no + 1, no-|-3^ ... no-|-l, no + L — 2<->no + l, no + L— 1 

I 

no -1-2, no + 3 

Notice that even though the other bound pairs far from 
the defects are also in resonance with the above states, 
the coupling with them occurs in second order in J, so 
they form a separated band. A more thorough discus- 
sion can be found in Ref.jT^, but the idea is the following. 
There are two possible channels for the second order cou- 
pling between distant bound pairs and the d-band states. 
A possibility is a transition from sites (no -I- 3, no -I- 4) 
to (no + I, no -f 4), but we can equivalently hop from 
(no -|- 3, no -I- 4) to the bound pair next to the defect 
(no-l-2, no-l-3) [two channels also exist for the excitations 
on sites (no — 3, no — 2)]. The amplitudes for these two 
transitions are equal in magnitude, but opposite in sign. 
Due to such quantum interference, or antiresonance, we 
have two uncoupled bands of very close energy. 

The antiresonance explains why the drop in the value 
of the maximum concurrence when A = d/J is not so 
drastic as one might have expected. Also, since there 



where fci = 1, 2 and ^2 = 1, 2, ...,L - N - A. 

Notice that the states above cannot be obtained when 
we have small chains with L < 8, 9 and the total number 
of excitations are N = L/2, L/2 + 1 respectively, as the 
available L — 6 sites are not enough for the remaining 

— 1 excitations. But, for these small chains and also 
for larger ones, there are other sorts of Bell-type states 
that can be obtained making use of the anisotropy. An 



are only four states preventing the creation of a perfect 
Bell state, the decay of Cmax becomes less perceptible in 
larger chains. 



4- Two excitations and d/J 

When A 3> d/J, the allowed transitions, shown below, 
now have two registers out of the core that could lead to 
a Bell-type state. 

no,no-|-2 no, no+3 ^ ... no, no + L—2 

I : 

no + l,no-|-3^ ... no-l-1, no + L — 2Mno-|-l, no-l-L-1 

This justifies why the maximum concurrence never goes 
back to 1, even for very large A. However, the effects 
caused by just these two isolated registers, 0(no,no -I- 2) 
and (/'(no -|- 1, no -l-L — 1), become less significant in larger 
chains. 

5. More than two excitations and A ^ d/J 

When more than two excitations are present, the anal- 
ysis of the maximum concurrence for small A's gets more 
complicated. A more interesting and general situation 
emerges when d/J. Here, we have again Bell-type 

states similar to the ones given by Eq. (|14|) . These states 
have one excitation hopping between the two defects and 
all the others bound together in a cluster in sites far from 
the defects. 

Only clusters with the same number of excitations are 
coupled and the transitions between them happen in a 
high order of perturbation theory in J. These transi- 
tions involve virtual steps where a dissociation occurs. 
Therefore the clusters move together as a whole and very 
slowly. We have found that states involving clusters at 
least two sites away from the defects have the largest 
concurrences. They are written as 



(27) 

I 

example is the state 1 / ^/2[(j){no , no + 2,no + 3,no + L — 
2, no-l-L — l)-|-(/)(no-l-l, no -1-2, no -1-3, no-l-L— 2, no-l-L— 1)] 
found when L = 10. We also note that, in the case of 
L = 12 and N = 5, the states with large concurrences 
obtained numerically did not correspond to what was ex- 
pected from H27|l . though they were also Bell states. 

As a result, the anisotropy coupling may be harmful 
when small, but very useful when sufficiently large, for 



no+L-N-l 

\^{n,mi,m2, ■■■,niN-i)) = ^ ^sin 

mi=na+4 

1 



7rA;2(mi — no — 3) 



{L-N-3) 

[(/■(no, mi,mi + 1, mi -f- A^ - 2) + (-l)'=i+V(no + l,mi,mi + 1, ...,mi+N - 2) 
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we can recover Bell-type states. 

In the limit of very large A we can use the forma- 
tion of clusters to create maximally entangled states that 
will remain as such for a long time. As an illustration, 
we take the case of L = 12, 4 excitations, d/J= 10, 
A = 50 and the two defects placed on sites 1 and 2. The 
most straightforward method to study the dynamics of 
the system consists of diagonalizing the 495 x 495 Hamil- 
tonian for the chosen parameters. We then assume 
that the initial state is the quantum register 0(1,6, 7,8) 
with one excitation on site 1 and the other three exci- 
tations bounded in a cluster on sites 6, 7, and 8, which 
is not an eigenstate of the Hamiltonian H. Its evolu- 
tion is obtained by writing it as a linear superposition of 
the eigenstates of H with their respective energies. The 
results are shown in Fig. [2 

The excitation on site 1 hops between the two defects 
with a period 2tt/ J The analytical reason for this is 
the following. The states with one excitation on one of 
the defects and three excitations bounded together form 
their own energy band around the value + d -\- JA. 
These states are not coupled with the other states of 
the system and can be treated separately, largely reduc- 
ing the Hamiltonian to be considered. In first order of 
perturbation theory, the Hamiltonian is still further re- 
duced to a simple 2x2 matrix written in the basis of 
the two registers (/)(1,6, 7, 8) and (/)(2, 6, 7, 8), which are 
coupled in first order. The off diagonal elements of this 
matrix are J/2. Trivially, the eigenstates and eigenvalues 
are the EPR states V-i = l/\/2[(/)(l, 6, 7, 8) ± 0(2, 6, 7, 8)] 
and E± = 64 + d + JA ± J/2, respectivelly. For the 
initial state considered, the probability to find the ex- 
citation on site 1 at time t is given by J0(i,6,7,8) (0 — 
{1 + cos[(£'+ — E^)t]}/2, while the probability to find it 
in 2 is P^(2,6,7,8)(0 = {1 - cos[(j;+ - E_)t]}/2. 

At each instant of time kir /2J , where k is an odd num- 
ber, a state of maximal entanglement between the defects 
is obtained, l/\/2[0(l, 6, 7, 8) ± 0(2, 6, 7, 8)] . After a very 
long time, the states 0(1,6,7,8) and 0(2,6,7,8) finally 
start mixing with the other six registers (see Eq. lj23l), 
where the cluster appears in a different position. The 
hybridization of these states in time is shown in Fig. |21 
The top panel corresponds to a short time and the two 
lower panels are obtained after much longer times have 
passed. 

Before allowing the system to evolve for a very long 
time, where the cluster delocalizes, states with concur- 
rence very close to 1 can be created. We let the system 
evolve just up to one of the initial instants where a Bell 
state emerges. At this moment, the two defects can be 
quickly detuned and such maximally entangled state be- 
tween the defects would be maintained for a long time 
until the cluster starts moving. 

It is clear that if the initial state was one of the EPR 
states l/\/2[0(l,6, 7,8) ± 0(2,6,7,8)], it would take a 
very long time to change, since this would require mov- 
ing the "heavy" cluster of three excitations together. 
The ability to keep the entanglement between the de- 



fects very large for a long time is a consequence of the 
large anisotropy of the XXZ model. In this way, this 
result strongly differs from the dynamics considered in 
Ref. There, such effects of the Ising interaction are 
not seen, because the dynamics of correlations is studied 
for an XY model. 
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FIG. 2: (color online) We have a chain with 12 sites, 4 excita- 
tions, and two defects placed on sites 1 and 2. The parameters 
of the system are d — lOJ and A = 50. We choose 0(1, 6, 7, 8) 
as the initial state and let the system evolve. The three pan- 
els show the probability for finding other quantum registers in 
time. For the short times of the top panel, 0(1, 6, 7, 8) (black 
solid line) mixes only with 0(2, 6, 7, 8) (red dotted Une). The 
instants of time where Bell-type states are created are shown 
with vertical long dashed lines. They appear at 7r/2J, 37r/2J, 
and 57r/2J. In the middle panel, maximally entangled states 
are not obtained anymore, as other quantum registers start 
mixing with the previous two. The probabilities for the regis- 
ters 0(1, 5, 6, 7) and 0(1, 7, 8, 9) almost coincide and are shown 
with dashed lines (blue and green, respectively). The same 
happens to the registers 0(2, 5, 6, 7) and 0(2, 7, 8, 9), which are 
indicated with long-dashed lines (violet and indigo, respec- 
tively). The lower panel is obtained for times even longer, 
where the quantum registers 0(1,8,9,10) and 0(2,8,9,10) 
barely start appearing. They are shown in the figure with 
dot-dashed (maroon) lines, but can hardly been seen. 



B. Next nearest neighbor defects 

The description of the dependence of the maximum 
concurrence on the anisotropy coupling when the defects 
are next nearest neighbors is very complex, especially 
when several excitations are present. Just like in the end 
of the previous subsection, we restrict our study to the 
case of A ^ d/J, where general remarks can be made. 

The states with the largest concurrences are similar 
to the states given by Ea. H27|) . They have one excita- 
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tion hopping between the two defects, though they are 
now separated, so the transition occurs in second order in 
J. The other excitations are bound together in clusters 
placed at least two sites away from the defects. Com- 
pared to the case of nearest neighbor defects, these re- 
maining A^— 1 excitations have now one site less available, 
there are L — 7 sites for them. 

When TV — 1 is larger than L — 7, the states with the 
largest concurrence have a different form. Sometimes 
these states can also have large concurrences, but this 
is not always the case, which explains the curves with 
low concurrence in the top and middle panels in Fig. |31 




FIG. 3: (color online) Maximum concurrence vs. anisotropy 
coupling in closed Heisenberg spin chains with defects located 
on sites 1 and 3. The energy difference between the defects 
and the other qubits is d = lOJ. The panels correspond to 
data for chains of different lengths: L = 8 (top), L = 10 
(middle) and L = 12 (bottom). The maximum number of 
excitations considered for each chain is L/2. The curves are 
chosen as follows: solid (black) curves for 2 excitations, dotted 
(red) for 3 excitations, dashed (blue) for 4 excitations, long 
dashed (green) for 5 excitations, and dot-dashed (violet) for 
6 excitations. 

The important result of this section is the verification 
that Bell-type states can again be obtained in the limit of 
large chains and A ^ d/J ^ 1, even though the defects 
are separated. 

An analysis for the time evolution of the system with 
a particular d-band quantum register taken as the initial 



state could again be made. It would be very similar to the 
one from the previous subsection, though a longer time 
is required for a maximally entangled state to appear. 
IV. CONCLUSION 

We have shown how the interplay between interaction 
and disorder may affect the entanglement between two 
defects in a chain of qubits with anisotropic coupling. 
When the difference in energy between the defects and 
the ordinary qubits d is much larger than the hopping 
integral J, several situations may be identified according 
to the value of the anisotropy coupling A, the number 
of excitations in the chain, where the defects are placed, 
and the chain size. 

The most general case refers to large chains and A ^ 
d/J » 1, where Bell-type states are obtained even for 
defects that are not nearest neighbors. These states have 
one excitation hopping between the two defects and all 
the others bound together in a cluster far from the de- 
fect sites. These clusters are a consequence of the large 
anisotropy. They move together and very slowly, allow- 
ing the maintenance of a large entanglement between the 
defects for a long time. 

In the case of nearest neighbor defects, we also verify 
that even when the anisotropy coupling is not very large, 
the decrease in the value of entanglement between the de- 
fects is never very abrupt and it becomes less significant 
in larger chains. Analytical results were obtained for two 
excitations when A = and A < d/J. For A < d/J 
the state with the largest concurrence is the one with the 
smallest energy in its band. 

We note that in this paper we have not developed any 
analysis in terms of quanturn phase transition as done 
in several previous works 

H Efl II3- By selecting two 
very large defects we guarantee that the entanglement 
between them is mostly kept very large. If only the XX 
part of the Hamiltonian was present, i. e. if A = 0, 
their concurrence would always have the maximum value 
1, as seen from Figs. ^ and O It is the effect of the 
extra Ising interaction that can sometimes decrease the 
entanglement. 

The present work has its clear relevance for quantum 
information and quantum computing, but it should also 
be of interest for condensed matter physics, where one 
wants to understand how localization may be affected by 
interaction and disorder. 
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